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STATEHEMT OF WORK l\~kO 

PRINCIPAL INVESTIGATOR: Dr. J. Ronald Burrows 

OBJECTIVES: 


The objective is to analyze disturbance fields caused by global external 
current systems in order to: 

a. Gain an improved understanding of the physical processes which 
control high latitude current systems. 

b. Increase t>s confidence level in the identification of internal 
(anomalous) field levels. 

APPROACH: 


The basic approach is to: 

a. Categorize the vector data by those physical parameters 
(gcomagneti c activi ty. Interplanetary conditions, etc.) important for 
investigation of external current systems. 

b. Map the disturbances for appropriate conditions. 

c. Model the currents which might cause the mapped disturbances. 

d. Correlate results with data from other sources. 

TASKS; 

a. Develop sub-sets of vector data,, from both polar regions, 
which are sorted by geomagnetic activity, interplanetary magnetic 
field conditions and solar illumination of the ionosphere. 

b. Investigate the distribution of vector residuals caused by 
external currents in the auroral oval and polar c.aps for various 
conditions. Some of the more immediate problems are: 

(1) Cusp-related disturbances during magnetically quiet 
condi tions. 

(2) Polar cap disturbances during northward interplanetary 
magnetic field (IMF) movement. 

(3) Interhcmisphcric asymmetries in polar cap disturbances 
as a function of IMF sector structure. 

{U) The dependence of davside disturbances in the northern 
hemisphere upon solar elevation angle. 
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c. Determine the spatial distribution of f ieid-ai i9ned current 
systoni$ in 'Mch data subset as a function of geomagnetic latitude and 
iocai magnetic time* 

d. Construct modeis of typical global ionosphere-magnetosphere 
current patterns for each subset, using ground based magnetometers 
and the derived distribution of field-aligned currents as data bases. 

e. Correlate the characteristics of inferred field-aligned 
currents with auroral arcs observed by the Defense Meterological 
Satellite Program (DHSP) satellites' optical and particle instrumentation 
on near-simultaneous high latitude transits of Magsat and the 

OMSP sate) 1 i te. 

The final report shall specifically describe, at a minimum, the results 
of each of the above tasks. The spec! fic methodologies and inter- 
pretation techniques (including programs) applied to address each 
shall be discussed in sufficient detail so that they can be cwnpletely 
understood. References to other publications are acceptable if their 
level of detail is adequate. If references as noted above are utilized, 
copies or excerpts from these references shall be provided along with 
the final report. The final report shall identify and discuss any 
significant results derived as a result of the Hagsat project. The 
report shall describe any new models constructed as a result of this 
project as well as describe how they were constructed, and how they 
differ from the existing model (s). 

1. Personnel 


Ors. T.J. Hughes, Margaret 0. Wilson, D.D. Wallis and the 
author have been involved in defining the investigative approach and 
developing analytical tools for the past year. Or. Wilson has been res- 
ponsible for data processing operations and software development; Dr. 

Hughes has been responsible for current modelling techniques and soft- 
ware; Dr. Wallis has been responsible for definition of coordinate sys- 
tems and related analytical work, development of* a data base of correla- 
tive magnetic rreasurements and establishing contact with ground-based 
magnetic and optical observations. 

2 . Techniques 

2.1 Data Processing 

Data processing work has been directed towards developing 
programs for the reduction and presentation of data in a variety of formats. 
Routines have also been developed to diagnose problems on the source 
(Chronicle) tapes and in the processing software. 

Data processing programs are designed to use CHROilFIN or 
CHPONINT tapes as input. Standard primary processing generates a WGEDT 
output tape which is used as input both for plotting routines and further 
analysis. This tape contains data records of 62.9 seconds duration in 
which the spacecraft position (in geocentric coordinates), scalar and vec- 
tor residuals and model field (resolved in NEV components) are computed 
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at second resolution. CHRONFIH vector data is filtered with a 17 point 
Bartlett filter and decimated to every 16th point before the residuals art 
computed by subtracting the ^FC fIGST (^/80 field model. For the X and Y 
residual components, a 'fluctuation parameter' is conpited which pn^vides 
a measure of the small scale structure present in the data prior to fil- 
tering and decimation. This parameter is connsuted about once per second 
(i.e., every 16 data points). 

2.2 Data Display 

A Hunger of formats for displaying the data have been develoiwd. 
Standard and variable format plots of three orthogonal components of the 
vector residuals versus universal time are most frequently ^nerated. 

The fluctuation paran«ters and scalar magnetometer residuals m^ also be 
selected for display in these formats. Standard plots are limited to 
poleward of 50® eccentric dipole latitude because of the principal interest 
of this investigation. However, the vriable format plots are unrestricted 
as to universal time or latitude. Polar plots of the residual vectors 
in the plane perpendicular to either the vertical direction or the model 
field direction may be produced with varying time intervals betMeen 
successive vectors. Figures 1 and 2 show a test pass plotted in universal 
time and polar coordinates. The positional coordinates and frame of the 
residuals are based on an eccentric dipole system. 

2.3 Coordinate Systems 

Studies of the correlation of external field residuals with 
geophysical parameters are strongly influenced by the coordinate system 
chosen for the analysis. Considerable investigation of different coordinate 
systems has been undertaken. Figures 3 and 4 show polar plots of the same 
test pass in two other coordinate systems. Figure 3 uses the HDZ usually 
used for ground-based measurements. In figure 4, the coordinate frame 
for the residuals uses the model field direction, U, and the vector to 
the eccentric dipole origin, Z'. Differences in the orientation of the 
residuals in different coordinate frames is evident, as expected. Work 
will continue to establish the best coordinate frame for data analysis. 

3. Accompl ishments 

3.1 Software Ctevelopment 

Data from both intermediate and fine attitude chronicle tapes 
have been selected for processing and display. Since the fine attitude 
tapes are significantly better for studies of vector data, it has been 
decided not to use significant computing resources to routinely process 
interirediate attitude tapes because this work will only be duplicated 
later with the fine tapes. The CHRONFIN test tapes provided have been 
used to develop and test software. A list of some of the programs and 
subroutines is attached as Appendix A. Although not complete, it indicates 
the types of capabilities that have been developed. 

3.2 Correlative Studies 

Various magnetic parameters of the magnetosphere and inter- 
planetary medium have been acquired for inclusion in a correlative data 
base. The AE, AL, AN indices are not yet available. As indicators of 



auroral electrojet activity. It is expected that they will be among the 
BKjst useful for correlative studies. 

Discussions with ground-based observers of visual aurora 
In Alaska and Scandinavia have Identified some observing periods when 
Magsat made close crossings. Cofl|)ar1sons of auroral morphology and currants 
In the auroral oval. Inferred from Higsat, will be Initiated when fine 
attitude data for these periods becomes available. 

Due to son« operational failures of DMSP satellites, the 
satellite coverage of optical aurora during Magsat's lifetime Is less 
than anticipated. A con^uter search for simultaneous observing periods 
has not yet found any occasions that are suitable for study. 

3.3 Modelling of Current Systems 

Soft\<are Is nw available for model 11 no of global Ionospheric/ 
magnetospheric current distributions poleward of 60° 5 »omagnet 1 c latitude. 
This forward model 11 no procedure Is based on Kisabeth's method of specify- 
ing current distributions In the Ionosphere In cells 0.5° In latitude and 
4*^ in longitude and confuting the resulting disturbance fields along the 
Magsat orbit due to these currents and the field-aligned currents which 
are necessary to maintain current continuity. Current distributions are 
adjusted to obtain optimal fits to the vector cunsjonents Of the residual 
field measured by i4agsat. An exarrole using data from the fine attitude 
test tape is shown in figures 5 and 6. Figure 5 shows the ionospheric 
current in each cell of the model. Figure 6 shows the disturbance resi- 
duals at 300 km altitude along the soacecraft track due to this current. 

4. Publications 


A paper (No. 3.4.6) entitled "Small scale current sheets in 
the piannetosphere near the auroral ionosphere" by J.R. Burrows, Margaret 
D. Wilson and D.D. Wallis was presented at the 23rd COSPAR Plenary meet- 
ing in Budapest in June, 1980. It was based on data from the CHRONINT 
test tape. The paper was not submitted for publication in the conference 
proceedings because fine attitude data was not yet available for that 
period by the submission deadline. 

5. Data Quality and Delivery 

A principal concern at present is the high frequency noise 
which is introduced into the fine attitude test tape data by the irethod 
that the Project Office has adopted for computing the fine attitude. 

Some n^thod for smoothing the individual attitude solutions needs to be 
developed before standard production of fine attitude tapes is started. 
Discussions with the Project Scientist are continuing and a satisfactory 
solution is likely to be reached. 

The early derivation and distribution of Magsat field models 
by the Project Office has been very uKJCh appreciated since they are the 
keystone for our data analysis and subsequent interpretation. 

6. Conclusions 


The Magsat data set appears to have remarkably high precision 
and quality. It should permit major advances to be made in modelling current 
distribution at high latitudes in the ionosphere and magnetosphere. 



APPENDIX A 


Data Processing and Diagnostic Pn)qraws 

MAGATT Prografln« to read the attitude flag data fr<m the Magsat Input taj^s and 
plot and/or list the data. 

MAGNUM Programne to dun^ data from Magsat chronicle Input tapes (CHRONINT, 
CHROUSC and CHRONFIN types). Under Input card control the program 
will function In a variety of nudes. 

MAGLST Programme to list the contents of a 11A(^DT output tape from orbit and 

decimated data records (the tape can be eltiier normal or extended type). 

MAS^AL Programme to plot data from the MAGEDT tapes created by TEftDT. Under 
Input control, either decimated or undecimated data can be plotted, 
residuals or vectors can be rotated, time and latitude ranges can be 
set and several sets of data can be plotted In one job. 

SHMITS Programme to confute the Schmidt normalized coefficients for one or more 
field models and produce, for each model, an object nudule containing 
these coefficients. 

SPLTST Programme to investigate the differences between the computed and 
spline fitted values of the nudel field conponents using different 
separations of the grid points for the spline fit. 

TEICDT Standard processing programme to create an edit tape, MAGEDT, from Magsat 
CHROtlINT or CHRONFIN input tapes. The 'MA(£DT tape will always contain 
records of filtered decimated data and optionally may contain records of 
unfiltered, undecimated data. 


Data Procissinq and Diagnostic Subroutines 


MAGFLD Assembler subroutine Mhich computes the model field components BR. BT« 

BP for a given spacecraft position. Called by TEHEOT, SPITST. 

MAGGET Subroutine to read data from Magsat chronicle tapes. Called by TEMEDT, 
f'*AGDUM, mGAH. 

FORHIT Assembler subroutine which formats the orbit data and the scalar/vector 
data for listing. Called by TEMEDT. ^GDUH. 

STDPRE Subroutine to prepare data for the plotting subroutine STDPLT of TEMEDT. 
North and south polar data (from 50° eccentric dipole latitude) are 
plotted versus Universal Time. The vector residuals are rotated to the 
fIDU frame. Scalar, vector and fluctuation data are plotted with an 
x-scale of 120 seconds/inch and a plot decinatlon length of 4. Called 
by TEf€DT. 

STDPLT Subroutine to produce the standard decimated data plot for TEMEDT. 

Data to be plotted have been transferred to the arrays by STDPRE. Called 
by TEMEDT. 

POLPLT Subroutine to plot the polar diagram for Magsat transpolar plots. 

Called by STDPLT, MAGPLT. 

MAGROT Subroutine to rotate Magsat NEV residuals to one of five other systems. 
Called by STDPRE. PREPAL. 

6E0CEN Subroutine to convert inertial coordinates to geocentric for 128 minutes 
of Magsat orbit data. Called by MAGGET. 

MAGSPL Subroutine to con^jute cubic spline interpolation coefficients for various 
functions y(x) using successive sets of 6 grid points and storing In 
the output array the coefficients at the third grid point in each set. 
Grid points are assumed to be 60.000 seconds apart. Called by TEFEDT, 
SPLTST. 

GETRTP Assembler subroutine computes polar coordinates R, T, P from cartesian 
coordinates X, Y, Z. Called by MAGGET, TEFEDT. 

ECOIP Subroutine transforms geocentric cartesian points to eccentric dipole 
cartesian coordinates. Called by TEfEpT. 

FILTM Subroutine to filter and decimate a Magsat data time series. Called by 
TEl-EDT. 

FLDTRC Subroutine to compute the position of a point in geocentric cartesian 

coordinates when traced down the magnetic field, specified by its direc- 
tion cosines, to a reference altitude above the earth's surface. Called 
by TEf'EDT. 


^TXYZ Subroutine to obtain geocentric cartesian coordinates of points having 
coordinates R, T, P where R Is the vector iwignitude either relative to 
the ^ocentrlc origin or to the earth's surface, T Is latitude In 
radians and P Is longitude In radians. Called by TEfCOT. 

PREPAL Subroutine to prepare data from tai^s for plotting, using either 

decimated data or extended (undeci mated) data. Called by HAi^Al. 

MAGPLT Subroutine to plot the flagsat residual data accuimila^d In arrays by 
PREPAL. Called by HASPAL. 

SYNTHS Subroutine to produce synthetic scalar data from the decimated X, Y, Z 
vector data. Called by PI^PAL. 

RECONS Subroutine to reconstruct measured vector data from model field and 
residuals stored in MAGEDT tapes, as required either for the rotation 
of these vectors to obtain fluctuation parameters In the rotated frame, 
or for the coroutation of synthetic scalar data from the corresponding 
X, Y, Z con?)onents of vector data. Called by PREPAL, WGPAL. 

FLUCTR Subroutine to cotmsute the X and Y fluctuation parameters from rotated, 

reconstructed, undecimated X and Y vector measurements. Called by PISPAL. 

POLVEC Subroutine to produce a polar plot of the residuals In the horizontal 
plane or plane perpendicular to the nradel field direction. The space- 
craft position is normally specified in eccentric dipole latitude and 
eccentric dipole magnetic local time. Called by MA6PAL. 


Current and Maqtutlc Perturbation Modelling Proqrawros 

CELSAM This code Is designed to read data from FT10» idrtilch has been ^nerated 
by DATPRP. The latter reads HAC£0T tape, and reformats the data for 
plotting using TJH's plot routines (to be used with no^l data}. 

The code converts all latitudes and longitudes to the comsponding 
cell positions (observer kernel positions), coos^utes the cell numbers, 
locates the data point closest to the cell observation (mint. The 
closest point Is the nrinliraim sum of squares of differences beteeen the 
real point and the cell point (longitude Is normalized to a distance 
via sin (col at). Output Is hard copy and FT11. 

EWFAOl Driver to examlr^ the Integration of the EH system uses Blot-Savart; 

Integrations are done by a Roidierg scheme applied to Intervals of dif* 
ferent sizes (Input). 

FILRES Programnm to confute the Impulse and/or the Heaviside response of recur* 
sive type filters. 

FILTER A general code for examining the frequency response of a nunter of fil- 
ter types (I.e., non- recursive, non-cascaded; non-recursIve, cascaded; 
recursive). 

1-DDEL4 This code generates kernels for use with the grid cell model of Kisabeth. 
This code is not generalized. It Is assumed that there are 90 cells 4 
degrees long (longitude) and 54 current, rings 0.5 degrees In latitude 
(58 to 85 degrees N) and that there are 77 observer positions (50 to 88 
degrees north In 0.5 degree steps. Maximum use of symmetry has been 
used, and only 46 kernels per current ring are actually conmuted and 
stored. Singly dimensioned variables are used throughout for efficiency. 
The magnetization term for the north-jouth current system perturbations 
is not calculated, but It is assumed that the Z-transform w1 11 look 
after this where appropriate. Output Is to 77 files, one for each obser- 
vation ring. Each file Is DDEFed with two sequence nun^rs (001 and 002) 
to accomnvjdate the NS (001) and EW (002) systems, so that In fact there 
are 154 separate files of kernels. 

nSDRIV Driver for north south current system 3 field perturbations using Kisa- 
beth 's equivalent charge method, and using Ron^erg Integration scheme 
exclusively. The code uses Kisabeth 's charge formulation to confute 
the B field due to a north south current system. This Involves Integra- 
tion along lines of charge whose c«ometry is that of dipole field lir»s 
placed at the ends (E/H) of the N-S current. The integrations along 
these lines is carried out in three steos, depending on the position 
of the observer. See "NSSUBZ'' and the oarameter "DEL" In that subroutine 
for the locations of the sub-intervals. The form of Rortjerg used in 
this code has been modified to carry out three Integrals slnwltaneously, 
for the sake of efficiency. However, each integral Is done only until 
convergence is reached. See R0fS3A and R0FB3B for details. 

NSGL2 The code uses Kisabeth' s charge fonrxilation to confute the B field due 
to a north south current system. This involves Integration along lines 
of charge whose geoTOtry is that of dioole field lines placed at the 
ends (E/H) of the N-S current. Gaussian (2 point) Integration is used 
for both integrals approx. 35 sub-intervals along the field line, and 
either 1 or 4 in latitude, depending upon the obs. position. 
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TFAVG pt'Qaramne to con^uti tht transfer function, nndulus of sanit, and phase 
of cascaded single averaging filters, and then the result of (tecliMtIon 
based on textbook by Blackman. 

ITILT This code Is desigmd to handle tilted ovals. In the sense that tiM 
current flows EH/NS In a system idiose pole Is offset towards midnight 
an amount "TILT" degrees. The oval Is specified by Its poleward border 
at PHI>90 degrees (dawn) (pollat) and the approximate eastern and ««stem 
boundaries (approximate because these change In the rotated system). 

The code defaults (ICUR»C) to constant current density everywhere with 
d(NS)*0.5* cur aaq)/m, and J(EH)a-1.0* cur amp/m. Otherwise subroutirw 
’CURRENT" must exist, and Is used to compute the current densities for 
each cell. Observer positions are rotated from the tilted oval system 
to that in which the oval Is a centered circle (the primed system). 

The closest positions for which there are computed kernels are then 
found, and these are reverse transformed to the tilted system for 
future reference, plotting, etc. Tne defined current densities are 
transformed to the primed system where the appropriate totals are com* 
puted for each cell. Finally, H', O', Z' are projected Into the obser- 
ver frame to give H, D, Z due to a tilted current system. 


Current and Magnetic Perturbation Model lino Subroutines 

CHGFUN Computes the integrands for the integrals used in Kisabeth's equivalent 
charge method for calculating the magnetic perturbations due to a north 
south current system. The integrand is computed as a function of obser- 
ver position, field line - ionosphere Intersection position and position 
along the field line. Mo integrals are carried out. Along the field 
line, resolution in colatitude is variable as specified below. The 
distance in kilometres from the observer to the position on the field 
line is also confuted. 

CLOSEST Subroutine to compute the spherical geocentric co-latitude angle des- 
cribing a point of a specified field line which Is closest to a given 
observer position. 

R01B3 This code computes three integrals siwltaneously. If Jl, J2, or J3 is 
equal to zero, the corresponding integral is not calculated. The code 
is ideal for use in calculating the three conxJonents of the magnetic 
field, where many of the factors and terms are the same for eaci* con^jo- 
nent. Vihen convergence in a component has been achieved, the appro- 
priate flan is set to zero so that further calculations on this convonent 
are not carried out. High precision and rapid convergence depend on 
the existence of all derivatives in the closed integration interval 
(A, B). Any singularities will cause a poor answer, and waste coit^uter 
time. The advantage of this method over Gauss-Leoendre is that the 
precision of the solution is easily estimated. For a description of 
the method, see A2ST Course in Numerical Analysis', Ralston, 1965. 
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